We are IntechOpen, 
the world’s leading publisher of 


Open Access books 
Built by scientists, for scientists 


5.300 130,000 150M 


ailable International authors and editors Downloads 


Our author among the 


154 TOP 1% 12.2% 


Countries delivered to most cited s Contributors from top 500 universities 


Selection of our books indexed in the Book Citation Index 
in Web of Science™ Core Collection (BKCI) 


Interested in publishing with us? 
Contact book.department@intechopen.com 


Numbers displayed above are based on latest data collected. 
F information visit www.intechopen.com 


ay 


Chapter 1 


New Insights and Horizons from the Linear Response 
Function in Conceptual DFT 


Paul Geerlings, Stijn Fias, Thijs Stuyver, Paul Ayers, 
Robert Balawender and Frank De Proft 


Additional information is available at the end of the chapter 


http://dx.doi.org/10.5772/intechopen.80280 


Abstract 


An overview is given of our recent work on the linear response function (LRF) x(r, r') and 
its congener, the softness kernel s(r, 1’), the second functional derivatives of the energy E 
and the grand potential Q with respect to the external potential at constant N and p, 
respectively. In a first section on new insights into the LRF in the context of conceptual 
DFT, the mathematical and physical properties of these kernels are scrutinized through 
the concavity of the E = E[N, v] and Q = Q [u, v| functionals in v(r) resulting, for example, 
in the negative semidefiniteness of 7. As an example of the analogy between the CDFT 
functionals and thermodynamic state functions, the analogy between the stability condi- 
tions of the macroscopic Gibbs free energy function and the concavity conditions for Q is 
established, yielding a relationship between the global and local softness and the softness 
kernel. The role of LRF and especially the softness kernel in Kohn’s nearsightedness of 
electronic matter (NEM) principle is highlighted. The first numerical results on the soft- 
ness kernel for molecules are reported and scrutinized for their nearsightedness, reconcil- 
ing the physicists’ NEM view and the chemists’ transferability paradigm. The extension of 
LRF in the context of spin polarized conceptual DFT is presented. Finally, two sections are 
devoted to ‘new horizons’ for the LRF. The role of LRF in (evaluating) alchemical deriva- 
tives is stressed, the latter playing a promising role in exploring the chemical compound 
space. Examples for the transmutation of Nz and the CC — BN substitution pattern in 2D 
and 3D carbocyclic systems illustrate the computational efficiency of the use of alchemical 
derivatives in exploring nearest neighbours in the chemical compound space. As a second 
perspective, the role of LRF in evaluating and interpreting molecular conductivity is 
described. Returning to its forerunner, Coulson’s atom-atom polarizability, it is shown 
how in conjugated n systems (and within certain approximations) a remarkable integral- 
integrand relationship between the atom-atom polarizability and the transmission proba- 
bility between the atoms/contacts exists, leading to similar trends in both properties. A 
simple selection rule for transmission probability in alternating hydrocarbons is derived 
based on the sign of the atom-atom polarizability. 
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1. Introduction 


A continuous challenge for theoretical and quantum chemists is to see if ‘classical’ chemical 
concepts describing bonding, structure and reactivity—the common language of all chemists 
—can still be retrieved from the nowadays extensive and complex computational results 
obtained at different levels of complexity with wave function or density functional theory. 
Conceptual density functional theory (CDFT) [1-6] has played an important role in this 
endeavour in the past decades. CDFT is a branch of DFT [7, 8] aiming to give precision to 
often well-known but sometimes vaguely defined chemical concepts (e.g. electronegativity, 
hardness and softness), affording their numerical evaluation, and to use them either as such 
or in the context of principles such as Sanderson’s electronegativity equalization principle [9] 
or Pearson’s hard and soft acids and bases principle [10]. ‘Chemical’ DFT or even ‘chemical 
reactivity’ DFT would have been a better name for the obvious reason that concepts are 
essential for all branches of DFT (especially the fundamentals) and that chemical reactivity is 
one of the main issues addressed in conceptual DFT. 


When looking at the basics of CDFT, the energy functional, E = E[N, v], stands out [8]. Why? It 
is the key ingredient to get (qualitative and quantitative) insight into the eagerness of an atom, 
or a molecule, to adapt itself to changes in the number of electrons, N, and/or the external 
potential, u(r), that is, the potential felt by the electrons due to the nuclei. These changes are 
essential in describing (the onset of) a chemical reaction, hence chemical reactivity. The readi- 
ness of a system to adapt itself to these new conditions is quantified through response functions, 
6” (0"E/ON") /dv(r1)...60(tm), which are the cornerstones of CDFT. Literature on these response 
functions is abundant, especially on the first-order responses (electronic chemical potential u 
[11] (n =1,m = 0) and the electron density p(r) (n = 0,m = 1), the cornerstone of DFT itself) 
and two of the second-order responses (chemical hardness ņ (n = 2,m = 0) and its inverse, the 
chemical softness S [12], and the electronic Fukui function f (r)(n = 1,m = 1) [13]). The most 


prominent of the third-order response functions [14] is the dual descriptor f Q) (r) [15], the 
N-derivative of the Fukui function (n = 1,m = 1). Remarkably, response functions diagonal 
in v(r) (n = 0, m = 2 — 3) were nearly absent in the CDFT literature until about 10 years ago 
(see [16] for an overview of this early work). The reasons are obvious; here we concentrate on 
its simplest member (n = 0, m = 2), the linear response function (LRF). 


x(r,r') = (E/Sv(r)óv(r))y (1) 


The calculation of this kernel turns out to be far from trivial, as is the representation of this 
quantity, a function of six Cartesian coordinates, and by extension its link to ‘chemical’ concepts. 
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Note that in the context of time-dependent DFT [17-19], the LRF has made its appearance 
many years ago as it was realized that the poles of its frequency-dependent form are nothing 
other than the electronic excitation energies. Thanks to Casida’s elegant matrix formalism [20], 
electronic transition frequencies, intensities and assignments are nowadays routinely performed, 
implemented as they are in standard quantum chemistry packages. However, this evolution was 
not accompanied by a parallel endeavour on the evaluation, representation and chemical inter- 
pretation of the frequency-independent or static LRF. 


In the past decade, the ALGC group, in collaboration with colleagues from different countries 
(Canada (Ayers), US (Yang), Spain (Sola), Poland (Balawender), etc.), set out a program aiming 
at the systematic evaluation, representation and interpretation of the LRF with the following 
results obtained until 2013, summarized in a review paper in Chemical Society Reviews [16] 
(no explicit reference to each of the individual constituting studies will be given). 


1. The LRF can now be routinely calculated at several levels of approximation for which the 
coupled perturbed Kohn-Sham perturbational approach turns out to be the most attractive 
approach, also permitting different levels of sophistication depending on the treatment of 
the exchange correlation potential (v,.) term in the perturbation equations. In its simplest 
form (neglecting the influence of the external potential variation on the Coulomb and 
exchange-correlation terms in the perturbational equations), the independent particle 
expression, already presented by Ayers and Parr [8, 21], is retrieved. 


2. The representation can be done via contour diagrams (fixing, e.g. r’) as demonstrated for 
atoms and molecules, or in the case of molecules, after condensation, via a simple atom- 
atom matrix, reminiscent of reporting the results of a population analysis. 


3. An abundance of chemical information was shown to be present in the LRF ranging from 
the shell structure of atoms, to inductive and mesomeric effects, electron (de)localization 
and (anti)aromaticity in molecules. 


In the present chapter, a synopsis is given of the progress made since then by the ALGC group 
in collaboration with other groups as witnessed by two of the authors (P.A and R. B.), both on 
fundamental and applied aspects, that is, on new insights into the properties of the LRF and on 
new areas where the LRF is at stake. In Section 2 the mathematical/physical properties of the 
LRF are revised together with those of its congener, the softness kernel s(r,1’), the latter 
playing a fundamental role in scrutinizing Kohn’s nearsightedness of electronic matter 
(NEM) principle. The extension of LRF in the context of spin polarized DFT is also addressed. 
In Section 3, we highlight the importance of the LRF in the emerging field of alchemical 
derivatives when exploring the chemical compound space. We illustrate the potential of 
alchemical derivatives in exploring the CC — BN isoelectronic substitution in 2D and 3D 
unsaturated carbocyclic molecules: benzene and the C¢o fullerenes. Finally, in Section 4 we 
show how the LRF (in fact its forerunner, Coulson’s atom-atom polarizability) can be used to 
predict/interpret the conductivity behaviour of unsaturated hydrocarbons, thus entering the 
vibrant field of molecular electronics. 
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2. Theoretical developments 


2.1. On the negative and positive semidefiniteness of the LRF and the softness kernel and 
thermodynamic analogies 


The properties of the LRF x(r, 1’) are intimately related to the concavity/convexity properties of 
E[N, v] that we addressed in recent years [22, 23]. Where E[N, v] is convex with respect to (w.r.t.) 
N [24], it is well established that E[N, v] is concave w.r.t. v(r) following the Jensen’s inequality. 


E[N, Av, + (1 —A)w] > AEN, vı] + (1—A)E[N, v2} A Ae[0,1] (2) 


(see for example Lieb [25], Eschrig [26], and Helgaker et al. [27]). In Figure 1 (after Helgaker 
[28]) we illustrate the physical interpretation of this concavity property. For a given potential 
vı, the associated ground state energy is given by the expectation value (1,|H[{vi]|),), with 
ground state wave function tp; the red point in the figure. Changing vı to v2, with constant wp, 
induces changes in energy linear in ôv(r) as a consequence of the term | p(r)v(r)dr in the DFT 
energy expression. Relaxing the wave function yields an energy lowering (blue arrow) until 
the energy obtained by applying the variational principle with v = v2 is reached (with associ- 
ated wave function w,). Consequently, the true energy (black line) will always be found below 
the tangent line ensuring concavity. 


E[v] 


(W1|A[v2]|V1) 


(W2|H [v2]|2) 


(Wi |A [vi] Y} 


Figure 1. Illustration of the concavity of the E = E[v] functional (after Helgaker [28]) (Reprinted by permission from 
Springer Nature, Copyright 2016 [22]). 
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A direct consequence is that the LRF is negative semidefinite. 
{| x(t, 1) O(r)O(r')drdr’ <0 (3) 


where O(r) is any continuous function [23]. This inequality shows up when considering the 
second-order variation of the energy 


dE?) = ¥, || (S°E/d0(r)50(r’)) .,60(1)d0(x' )drdr' (4) 


When adopting the ôv(r) = Vod(r — r”) choice for dv(r), where Vo is a constant one gets 


dE) = V3 x(t", x’) <0 (5) 


showing that the diagonal elements of the linear response function x(r, 1’) should be negative 
or zero. This result links the concavity of the E[N, v] functional to the diagonal elements of the 
linear response function defined on R? x R°. This negativity was retrieved in all our numerical 
results (non-integrated and condensed) (see for example [29, 30]). Its physical interpretation is 
straightforward through the definition of x(r,1’) as (ôp(r)/ðv(r'))y: if the potential at a given 
point r is increased (made less negative), this electron-unfavourable situation will lead to 
electron depletion at that point, yielding a negative y(r, r). 


In Section 4, we will point out that Coulson’s atom-atom polarizability rag [31] can be 
considered as a Htickel-theory forerunner of the LRF and exploit its properties in discussing 
molecular conductivity. Defined as 


TTAB = 0q,/0aB = OE /Oa40ag = TIBA (6) 


the analogy emerges as q; is the 7t-electron charge on atom A, whereas (the change in) the 
Coulomb integral ag is equivalent to a change in the external potential. The case A = B was 
explored by Coulson using his complex integral formalism for Hiickel’s n energy and its 
derivatives, proving that 74,4 can be written as 


+20 


naa = (1/7) | (Aaa (iy)/A(iy))?dy (7) 


—oo 


Here, A represents the Hiickel secular determinant (or characteristic function) and A44 is its 
counterpart with row A and column A deleted. It turns out that A44 (iy)/A(iy) is imaginary and 
therefore 7744 is negative (see also Section 4). 


When introducing the E = E[|N,v] functional, it was stressed that in the LRF the second 
functional derivative is taken at constant N. In the remaining part of this chapter, it turns out 
that no less important role in CDFT is played by the softness kernel, which is the second 
functional derivative with respect to v(r) of the Q = Q|u,v] functional (the grand potential) 
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at constant u. Here, we switch from the canonical ensemble to the grand canonical ensemble 
[32] connecting two ways of specifying the same physics via functionals of each time two 
variables differing in one pair that is connected through a Legendre transformation (see for 
example [33]). 


Ofu, o] = E[N [u, 0]; 0] — uN [u,v] (8) 


where u and N are conjugate variables which are related by identities N [u,v] =N and 
LN, v] = u. Its apparent analogy with classical thermodynamics will be addressed later. 


In [23], we pointed out that Q[u,v] just as E[N,v] is concave w.r.t. v(r), implying that its 
second functional derivative at constant u, the softness kernel (note the negative sign in the 
definition) [34]. 

s(r, r’) = — (6O /õv(r)ôv(r’)) , (9) 
is positive semidefinite. This property fits the well-known Berkowitz-Parr relationship [33] 
linking x(r, r’) and s(r, r’) (vide infra). 


s(r, r) = =x, r) +f )/n (10) 


where f (r) is the Fukui function and ņ is the hardness. Indeed, 


[|se nomoa = || xer) +da 
a1) 


_ {| x(t, DOMO )drdt + 1/1 | (f(2)6(x))2dr>0 


since the hardness is nonnegative and x(r, 1’) was shown to be negative semidefinite. 


The properties of x(r,r') and s(r,r’) incited us to reconsider the analogy between the DFT 
functionals E and Q [32] on the one hand and the macroscopic thermodynamic state functions 
U=U[S,V], F= F|[T, V], H = H|S,P] and G = G|T,P] on the other hand (internal energy, 
Helmholtz free energy, enthalpy and Gibbs free energy written as functions of volume (V), 
temperature (T), pressure (P) and entropy (S)). Parr and Nalewajski extended the notion of 
intensive and extensive variables (T,P) and (S, V), respectively, in thermodynamics to the 
variables in DFT functionals by classifying external variables as properties additive with 
respect to any partitioning of the electron density p(r) = p,4(t) + pg(r) [32]. In this way, p(r) 
and N are clearly extensive, and u and v(r) are intensive. The analogy between G[T, P] and 
Q|u,0| can now be stressed: both the state function G and the DFT functional Q contain two 
intensive variables. This situation leads to a remarkable property when formulating a DFT 
analogue of the stability analysis in macroscopic thermodynamics [33, 35, 36]. Concavity for 
G[T, P] in all directions then implies that 


d’G = Grr(AT)* + 2GrpATAP + Gpp(APY <0, (12) 
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where Grr, Grp and Gpp are the second derivatives of Gibbs free energy, (Gxy = oG/ OXOY). 
This stability condition implies the negative semidefiniteness of the Hessian matrix. It yields 


Gpp = (OV /OP)- = —KrV <0 (13) 
Grr = —(05/0T)»p = —Cp/T<0 (14) 


that is, the isothermal compressibility xr and the heat capacity at constant pressure Cp should 
be positive, and the condition that 


(GrrGpp — Gpr) = (krVCp/T — œ? V?) 20 & KrCp/T2a°V (15) 


where a = (OV/OT),/V is the coefficient of thermal expansion. Another classical example of 
such stability analysis is the entropy written as S[U, V] (at constant number of particles) where 
both variables are now extensive [33]. 


Let us now consider the analogy with Q |u, v]. It is well known that 
(A/d), = —S[u, v], (16) 


where S is the global softness [3]. As the rh.s. of (16) is negative, concavity for Q{u,v] in u 
shows up. As discussed above, Q is also concave in v(r), leading to the positive semidefi- 
niteness of s(r, 1’) (these expressions being the counterparts of (13) and (14)). The condition for 
concavity in all directions leads after some algebra (see [23]) to the condition 


(@2/),Au +i (09/00), /dp) Avtr)ar) +i (a /12),(P°Q/60(r)50(")) , 
= (asa /5v(t)), /2u) (aa /dv()), /du) \Av(r)drAv(r)dr 20 (17) 


and finally to 


{| (s(x, x)S — s(r)s(r’))Av(r)drAv(r' )dr’ = 0, (18) 


the analogue of (15) where the local softness 
s(t) = (8(60/d0()),/9u) = (0/4), (19) 


has been introduced [37]. Taking again for Av(r) and Av(r') Dirac delta functions 6(r — r”) and 
6(r — r”), one obtains the condition 
s(r,r)S>(s(r))? 20. (20) 


This inequality shows that the diagonal elements s(r, r) should be positive, as could be inferred 
from the concavity of Q[u,v] but, more importantly, they impose a restriction on the relative 
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values of the three softness descriptors s(r,1’),S and s(r) in analogy to the thermodynamic 
relationship between kr, Cp and a at given T and P (V = VT, P]). This result is compatible 
with the aforementioned Berkowitz-Parr relationship. Indeed, starting from their expression 
for r = r (see (10), the relations s(r) = f(r)/n) 


s(r,r) = —x (r,r) + s(r)s(r)/S (21) 


and knowing that x(r, r) <0 (vide supra) one obtains 


s(r, r) — s(r)s(r)/S 20, (22) 


retrieving our conclusions above. 


2.2. Kohn’s nearsightedness of electronic matter revisited 


We now report on our recent explorations [38] on Kohn’s NEM principle. Kohn introduced the 
NEM concept in 1996 [39] and elaborated on it in 2005 with Prodan et al. [40]. In his own 
words, it can be viewed as ‘underlying such important ideas as Pauling’s chemical bond, 
transferability, and Yang’s computational principle of divide and conquer’ [40]. Certainly in 
view of the two former issues, this principle, formulated by a physicist, touches the very heart 
of chemistry and so, in our opinion, it was tempting to look at it with a chemist’s eye. Why 
however is this issue addressed in this chapter; in other words, what is the link between the 
LRF and nearsightedness? 


The quintessence of the NEM principle is as follows (Figure 2): consider a (many) electron 
system characterized by an electron density function p(r) with a given electronic chemical 
potential u. Now, concentrate on a point ro and perturb the system in its external potential v(r) 
at point 1’, outside a sphere with radius R around rp at constant electronic potential u. Then, the 
NEM principle states that the absolute value of the density change at ro, |Ap(ro)|, will be lower 
than a finite maximum value Ap, which depends on rp and R, whatever the magnitude of the 
perturbation. As stated by Kohn, anthropomorphically the particle density p(r)cannot ‘see’ any 


perturbation v(r) beyond the distance R [ro, Ap] within an accuracy Ap, that is, the density shows 
nearsightedness. Kohn offered evidence that in the case of 1D ‘gapped’ systems (i.e. with 
hardness 1 larger than zero) the decay of Ap as a function of R (i.e. upon increasing |r — ro|) is 
exponential and that for gapless systems it follows a power law. This suggests that in the 
molecular world, where 1) is observed to be always positive, the electron density should only be 
sensitive to nearby changes in the external potential. In [38], we provided the first numerical 
confirmation of this nearsightedness principle for real, 3D, molecules. 


Again, why address this issue in this LRF chapter? Going back to Kohn’s formulation quintes- 
sentially a change in density at a given point Ap(ro) in response to a change in external 
potential Av(r’) at different points is analysed. These are the typical ingredients of the LRF 
(change in v produces a change in p), and in this case the process is considered at constant 
electronic chemical potential u, in order words (dp(r)/dv(1’)),, is the key quantity. This is 
nothing else than the softness kernel s(r, r’) with a minus sign in front. Indeed, 
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vir’) 


R/ 


1 \ =cte 
P(to) + 5p; (1o) s | i S 
u= cte 
vír’) + bv2(r’) 
Sv,(r’) > dv, (r’) 


) 


v(r') + 6v,(r’) 


i R 
gr 
P(ty) + bP 2(ro) D 
i R 8 
Fie. 
ols) + Spalt) ka pac v(e’) + vlr’) 
8v,(r’) > dv,(r’) 


——* palo) = 6p 2(ro) = AP (ro R) 


Figure 2. Pictorial representation of the nearsightedness of electronic matter principle: when ôv(r')2ôv:(r'), with r 
outside a sphere with radius R around ro, 6p9(ro) no longer increases—6pp(to) = Ap (ro, R) no matter how large dv(r’). 
The intensity of the red and blue regions represents the magnitude of v(1’) and p(r), respectively. 


(Sp(r)/dv(r')),, = —s(t, r) = (6° Q/d0(x)50(7')) , (23) 


The Berkowitz-Parr relationship [34] can then be written as 
(5p(r)/S0(t'))y = (p(t) /50(r')),, — (dp(r)/du), (544/50(2')) y (24) 


an equation transforming conditions of constant N into constant u, for taking the functional 
derivative of p(r) w.r.t. dv(r) in analogy with this type of equation for partial derivatives in 
macroscopic thermodynamics [33]. As (dp(r)/ðu) , equals the local softness s(r) and (5/dv(r’)) , 
is an alternative way, as compared to (dp(r)/ON),, to write the Fukui function f(r), one retrieves 
from (28) 


slr, r) = -xlr r) + f@yf(r)/n (25) 


which will be the key equation in this section. As analytical methods are available to evaluate 
x(r,r'), f(r) and 7 on equal footing [41], s(r,r') can be evaluated, and its difference with x (r, 1’) 
analysed, in order to scrutinize nearsightedness. 
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(1,3) (1,5) 
“1.4 l 
a4 = linear response — softness ft 
As — softness f“ — softness f° 
-1.4 
“1.5 (1,1) 
-4.6 


Figure 3. Atom condensed linear response function and softness kernel of 1,3,5-hexatriene. The curves of the softness 
kernels using f+ and f— are overlapping [38]. 


In Figure 3, we depict the atom condensed linear response function and the softness kernel, the 
matrices X4g and sag, of 1,3,5-hexatriene. Condensation was performed by integrating the 
kernels over the domains of atoms A and B, represented by V4 and V4, that is, 


SAB = | | st r )drdr’ (26) 


Va Vg 


From our previous work on polyenes [42], the X4g matrix elements are known to show an 
alternating behaviour with maxima on mesomeric active atoms (C2,C4 and C6) and minima for 
mesomeric passive atoms (C3 and Cs), with the change in v(r) taking place at C1.The picture 
illustrates that the softness kernel is more nearsighted than the LRF (all its values are lower) 
and that sı,5 and s1,6 are very close to zero; the effect of the perturbation has died off completely, 
confirming the nearsightedness of the softness kernel. This effect can be traced back to a cancel- 
lation of the LRF, which is non-nearsighted, by the second term in Eq. (25), and accounts for the 
density changes induced by charge transfer from the electron reservoir to keep the chemical 
potential constant. Note that this condition, at first sight somewhat strange, is an often more 
realistic perspective when considering for example the reactivity of molecules in solution where 
the chemical potential is fixed by the solvent allowing (partial) charge transfer to or from the 
molecule while keeping its chemical potential constant [43, 44]. 


As a second example, we show in Figure 4 the change in density of the 1,3,5-heptatrienyl 
cation when the C atom of one of the terminal CH? atoms was alchemically replaced by an N 
atom. The corresponding density difference was evaluated through the alchemical derivatives 
approach (Section 3) where the carbon atom of the CH2 group was annihilated and replaced by 
a nitrogen atom at constant geometry and constant number of electrons. It is clear that the 
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wis’ 


Figure 4. Alchemical change in density using the linear response (top) and softness kernel (bottom) for the heptatrienyl 
cation to 1,3,5-hexatriene-1-amine [38]. 


effect of functionalization dies off much more quickly under constant u conditions. In fact, the 
third carbon atom is hardly affected in this—it should be stressed —unsaturated system, offering 
the possibility for mesomerism. In the LRF, the decay is slow and effects are still relatively important 
five bonds away from the perturbation. 


Further case studies on ‘3D’ systems (e.g. alchemically changing methylcubane to fluorocubane 
and on functionalized neopentane) yield similar results. All together, the results on the near- 
sightedness of the softness kernel found for all systems discussed in [38] are the first and a firm 
numerical confirmation of Kohn’s NEM principle in the molecular world. To put it in chemical 
terms, these findings provide computational evidence for the transferability of functional 
groups: molecular systems can be divided into locally interacting subgroups retaining a similar 
functionality and reactivity that can only be influenced by changes in the direct environment of 
the functional group. Thus, the physicist’s NEM principle and the chemist’s transferability 
principle [45]—at the heart of, for example, the whole of organic chemistry [46]—are reconciled. 


2.3. Extension of x(r, 1’) in the context of spin polarized CDFT 


As a natural extension of CDFT to the case where spin polarization is included [47-49], spin 
polarized conceptual DFT was introduced by Galvan et al. [50, 51] (see also Ghanty and Ghosh 
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[52] and for a review see [53]). In the so-called (Na, N B) representation, the response of the 
electronic energy to perturbations in the number of a electrons, Na, the number of f electrons, 
Ng, the external potential acting on the a electrons, v,(r), and the external potential acting on 
the £ electrons, vg(r), is studied. In this representation, the E = E[N, v] functional from Section 
2.1 is generalized to the E = E|Na, Ng, va, vg] functional. In the second, equivalent, (N, Ns) 
representation, in fact the one introduced by Galvan, the E = E[N, v] functional is generalized 
to E = E[N, Ns, v, vs] or E = E[N, Ns, v, B] where Ns denotes the electron spin number defined 
as the difference between the number of a and f electrons, Ns = Na — Ng, whereas v, is given 
by vs(r) = (va(r) — vg (4) /2. 


v(r) is equal to(va(r) + vg(r)) /2 and B represents an external magnetic field. If the magnetic 
field is static and uniform along the z axis, one has vq(r) = v(t) + pBz and vg(r) = v(t) — UpBz, 
v(r) is the usual external spin-free potential (as used in CDFT without spin polarization), 
whereas vs(r) is related to the magnetic field B(r). In this context, and sticking to the (N, Ns) 
representation, three linear response functions can now be defined: 


Xynnlr, r) = (© E/O0(4)50(2')) y Ne vg (27) 
Xsg (0,1) = (6°E/Sus(r)50s(t’)) y neo (28) 


Ans(t, r’) = (0° E/50(1)50s(t')) y, y, = (SP (8)/505 (1) 1, Nso 


. (29) 
= (dps(t')/50(t)) yy neo = (8 E/505(¥')50(t)) x, y, = Xs (r,r) 
where py(t) = Pa(t) + pg(t) and ps(r) = pa (£) — pg(r) are the total and spin densities, respec- 
tively. ynn (t, 1’) is the analogue of the spin-independent CDFT expression for the LRF y(r, 1’) 
(Eq. (1)). 


Space limitations prevent us to go in detail on the results reported in [47, 48]. We only depict in 
Figure 5 the SPCDFT analogue of the contour plots for y(r,1’) for closed shell atoms as 
discussed in [47] and the Chem Soc Rev paper [16]. In Figure 5, we plot the LRF, in the 
(Na, Ng) representation, for the ground state of Li ((1s)* (2s)'). The structure of x,,(1,1’) and 
Xgp(x, r’) is similar to the x(r,1’) plots for closed shell atoms: a negative diagonal part (cfr [47] 
and Section 2.1) and alternating positive and negative parts for r or r’ = constant in order to 
obey the trivial equation (cf. [47]) 


[xede = | (Spy (a)/50()) at = (8N/60())y =0 (30) 


Xaa extends further away from the nucleus than xgg in line with the extra a electron with the 
higher principal quantum number (n = 2) extending farther away from the nucleus than 
the n= 1 electron. 7,3 looks similar to the x, plot for He (insert) but contracted more to the 
origin (due to higher nuclear charge). The Xag and Xag plots show some evident symmetry, 


Xap(t, T) = Xgq(t’, r), and show positive regions along the diagonal. A positive perturbation in 
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Figure 5. Contour plots of the radial distribution function of the spin polarized linear response function of Lithium in the [Na, 


Nf] representation. r is represented on the horizontal axis, r’ on the vertical axis [(a) Lithium Xaa, (b) Lithium xag, (c) Lithium 
Xaa» (d) Lithium xgg. In the insert, the x, plot for He (see text) (Reprinted from [48] with the permission of AIP Publishing)]. 
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the a external potential ôvą(r) will cause a depletion of electrons in the vicinity of the pertur- 
bation, the 6 electrons are not affected directly. However, the depletion in a electrons will 
influence the Coulomb potential and due to the lower electron-electron repulsion, an accumu- 
lation in f electrons in the region considered will occur resulting in a positive diagonal xq (r, r) 
value. The concentration of the Xag and Xag isocontours along the r and r’ axes can be 
interpreted when referring to the Xaa and Xgg plots: perturbing vg(r) at a distance r’ larger than 
3 a.u. clearly has no effect on the p density, the Coulomb potential, the overall density and 
consequently on the a density. This results in (6p,(r)/dvg(t")) = Xag(t, 1’) having zero values 
for r larger than 3 a.u. On the other hand, perturbing vg(r) close to the nucleus induces a 
change in the 6 density, with repercussion on the Coulomb potential and so on the a density 
farther away from the nucleus (on the r axis) even in regions where the $ density is no longer 
affected. All these features account for the ‘partial plane filling’ of the Xag and Xag plots with a 


‘demarcation’ line at 3 a.u. 


To close this section, we mention that once x(r, 1’) (or its counterparts in SPCDFT) is known, a 
local version of the polarizability tensor components, say &yy namely ay,(r) can be obtained by 
straightforward integration: 


y(t) =- [koxt r')y(r)dr' 81) 


An example is given in Figure 6 [48] where for the atoms Li through Ne the trend of the 
spherically averaged a(r) (3 (axx(r) + &yy(r) + azz(r))) is given. From 2 a.u. on, the trends in 
a(r) for Li up to Ne parallel the global polarizability, known to decrease along a period of the 
periodic table. At lower distances (preceding the valence region), inversions with even negative 
values occur. The present results are evidently important when, for example, disentangling 


reaction mechanisms where local polarizabilities, that is, in certain regions of the reagents, are 


Figure 6. Plot of the local polarizability a(r) of the atoms Li through Ne via CPKS (see text) (Reprinted by permission of 
the publisher (Taylor and Francis Ltd.) [48]). 
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at stake and not the overall polarizability. The parallel between the relationship between local 
and global softness (Section 2.1) is obvious. 


3. The role of the LRF in alchemical derivatives and exploring chemical 
compound space 


3.1. Context 


The LRF has recently been exploited when investigating Chemical Compound Space [54-56]. 
Chemists are continuously exploring chemical compound space (CCS) [57, 58], the space popu- 
lated by all imaginable chemicals with natural nuclear charges and realistic interatomic distances 
for which chemical interactions exist. Navigating through this space is costly, obviously for 
synthetic-experimental chemists and also for theoretical and computational chemists who might 
and should be guides for indicating relevant domains in CCS to their experimental colleagues. 
Doing even a simple single-point SCF calculation at every imaginable point leads to prohibitively 
large computing times (not to speak about bookkeeping aspects and manipulation of the com- 
puted data). A very promising ansatz was initiated by Von Lilienfeld et al. [59-63] in his alchemical 
coupling approach where two, isoelectronic molecules in CCS are coupled ‘alchemically’ through 
the interpolation of their external potentials (see also the work by Yang and co-workers on 
designing molecules by optimizing potentials [64]). At the heart of this ansatz are the alchemical 
derivatives, partial derivatives of the energy w.r.t. one or more nuclear charges at constant number 
of electrons and geometry. The simplest members of this new family of response functions are: 


the alchemical potential 


Hal[N, Z, R] = (OW/dZ4) y r = (E/dZA)y r + (OVnn/OZa) yen = UAT + (32) 


and the alchemical hardness 


niglN, Z, R] = (0 W/0Z40Zz) yp = (WE/0Z10Z5) y r + (O° Vnn/dZ40ZB) y p = Ng + Ag = Ry, (33) 


where W(N, Z, R] = E[N, v|Z, R]] + Vnn[Z, R] is the total energy of a system,R = (Ry, Rz,...) 
denotes constant geometry and Z = (Z4, Zg, ..., denotes the nuclear charges vector. The anal- 
ogy with the electronic chemical potential and hardness (see Section 1) is striking. Some years 
ago, one of the present authors (R.B.) presented a strategy to calculate these derivatives for any 
atom or atom-atom combination analytically and to use them, starting from a single SCF 
calculation on the parent or reference molecule, to explore the CCS of first neighbours, imply- 
ing changes of nuclear charges of +1 or —1 [65]. In this way, simple arithmetic can be used, 
starting from a Taylor expansion 


W([dZ] = WIN, Z + dZ, R] — WIN, Z, R] = DHA dZ, + 12> DMs ndZadZpz (384) 


instead of a new SCF calculation for each transmutant. 
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The position of the alchemical derivatives in CDFT was already mentioned: they are response 
functions, now related to a particular charge in external potential, namely the charge in one or 
more nuclear charges. As the second derivatives are ‘diagonal’ in these particular external 
potential changes, a direct link with the LRF can be expected. Using the chain rule, one easily 
writes 


nis = (WEIN, piZ, R]/0Z40Z)y p= if (öp) /5v(r))y Olr) /8Z4) y rlr) /dZs) y parade’ 
(35) 
= | [xe r’)(1/|r — Ra|)(1/|r’ — Re|)drdr’ 


indicating that the alchemical hardness is obtained by integration of the LRF after multiplica- 
tion by 1/|r — R4| and 1/|r’ — Rg]. The basic relationship of the alchemical derivatives and the 
LRF shows how, again, the LRF makes its appearance in sometimes unexpected areas (see also 
Section 4). In the case of the first derivative, the (6E/6v(r))), factor in the integrand simplifies to 
p(r)(see Section 1), yielding. 


uy = | p(2)/fe— Raldr (36) 


the electronic potential at the nucleus, well known as the electronic part of the molecular 
electrostatic potential [66]. 


3.2. Applications 


As a very simple example, we consider the transmutation of the nitrogen molecule. Five chem- 
ically relevant mutants can be generated (see Figure 7) as nearest neighbours in CCS (AZ = +1) 
and at constant number of electrons: CO, NOF, OF, CN” and oc The differences between 
‘vertical’ (i.e. exact, via two SCF calculations) and alchemical transmutation energies were 
evaluated at the B3LYP/cc-pCVTZ level (note the inclusion of additional tight functions ‘C in 


CO NO oC, 


w 3 r 
oy CN e N <& ON 
eg 
E NC* OC 
Oy, 


Figure 7. Transmutation of the nitrogen molecule to its nearest neighbours in chemical compound space (Reprinted with 
permission from [54]. Copyright (2017) American Chemical Society). 
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order to recover the (changes in) core-core and core-valence correlation). The mean absolute 
error was only 0.034 a.u., the N — COtransmutation being particularly successful with a 
difference in energy of only 0.004 a.u., that is, 2 kcal mol ~+, The transmutation into a neutral 
system results in the cancellation of the odd terms in the Taylor expansion. Of course, we do 
not claim ‘chemical accuracy’, yet the ordering of the energy of all transmutants came out 
correctly, indicating that the alchemical procedure (even when stopped at second order) is a 
simple, straightforward road to explore CCS for neighbouring structures. 


Similar conclusions could be drawn for transmutation of benzene, for example, by the substi- 
tution of CC units by their isoelectronic BN units. The replacement of a CC unit in an aromatic 
molecule by an isoelectronic unit BN has been shown to impart important, interesting elec- 
tronic, photophysical and chemical properties, often distinct from the parent hydrocarbon [67]. 
An in-depth study of all azoborines (Figure 8) C6H6 — Co-2,H6(BN),, (n = 1, 2,3) turned out 
to reproduce correctly the stability of all possible isomers for a given n value (3,11,3 for 
n=1,2,3), which is of importance for applications in graphene chemistry where the 
(CC), — (BN), substitution is a topic that has received great interest in recent years [68]. 


As a computational ‘tour de force’ , and passing from ‘2D’ benzene to ‘3D’ fullerenes, we 
recently explored the alchemical approach to study the complete (CC), — (BN), substitution 
pattern of Ceo, all the way down to (BN)39. Ceo —> Coo-2n(BN),,, 1 = 1, 2, ...,30, predicting and 
interpreting via ‘alchemical rules’ the most stable isomers for each value of n. This study is 
based on a single SCF calculation on Ceo and its alchemical derivatives up to second order, 
enabling each possible transmutation energy to be evaluated by simple arithmetic (the diago- 
nal elements of the alchemical hardness matrix are equal) 


AW, [N = (N1, ..., Nn), B = (B1, ..., Bn)] = nnii + sy Ga T nn) = > > Nay (37) 


i=l j>i i=1 j=1 


where AW, is the transmutation energy from Ceo to Ceo-27 (BN), and N/B is a vector of the 
carbon atoms replaced by the nitrogen/boron atoms. 


2D 3D 


Figure 8. CC-BN Substitutions in 2D and 3D unsaturated carbocyclic systems (number of isomers for the 2D case in 
parentheses). 
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In (37), it is seen that the linear term drops as u; is unique by symmetry for all atoms in Ceo and 
because Ł AZ ; = 0 for any transmutation. The alchemical hardness matrix n thus completely 
determines the substitution energy and pattern. To summarize the results (for an in-depth 
discussion see [55]), the study reveals that the correct sequence of stabilization energies for 
each n is retrieved by adopting an approach in which for each n value the problem is looked 
upon without prejudice of the n — 1 result (called the ‘simultaneous’ approach). The other, 
simpler method (we called it the ‘successive’ approach) was shown to fail already after n = 13, 
be it that at some n values identical values were obtained with both approaches, for example, 
for the Belt structure (n = 20). Needless to say that even the ‘successive’ approach might 
already be prohibitively demanding for standard ab initio or DFT calculations, let it be for the 
simultaneous method. The sequence of substitutions could be interpreted in terms of a number 
of ‘alchemical’ rules of which the two most important are when (referring to earlier work by 
Kar et al. [69]): (1°) hexagon-hexagon junctions are preferably substituted, in a way that 
minimizes the homonuclear (BB,NN) bonds and (2°) the higher stability is created by maxi- 
mizing the number of filled hexagons. For the subsequent, more intricate, rules we refer to [55]. 


Shortage of space prevents us to comment on our recent results on the evaluation of isolated atom 
alchemical derivatives up to third order with different techniques, from numerical differentiation 
(not discussed hitherto in this chapter), via the coupled perturbed Kohn-Sham approach as 
discussed before, to the March and Parr combined 1/Z and N~'/°E[N, Z] expansion ansatz [70], 
permitting a walk through the periodic table on the road to scrutinize periodicity effects [56]. 


4. The role of the LRF (or its forerunner, the atom-atom polarizability) in 
evaluating and interpreting molecular electronics 


4.1. Context 


In this section the role of the LRF in molecular electronics is highlighted [72] has been a vibrant 
area of research in recent years. An ever-increasing number of papers (both experimental and 
theoretical) studied the transport properties of typically organic molecules containing Tt-conju- 
gated systems and considering possible applications for incorporation in molecular electronic 
devices (MED) [73]. Most of these theoretical studies have been performed at a high level of theory, 
but this type of calculation does not always lead to simple insights into why some molecules will 
conduct and which will insulate, and how the positions of the contacts influence this behaviour. 


We therefore adopted a simple ansatz based on Ernzerhof’s source and sink potential (for 
details see [74, 75]) in Fowler’s tight-binding Hückel approach [76]. One thereby considers 
only the r electrons of the molecule and cuts the resonance effects between the contact and the 
molecule after the molecule’s nearest neighbour in the contact. In the so-called weak interac- 
tion limit (see details in [71]), the transmission probability at the Fermi level T(0), which is 
directly proportional to the conductance, can then be written as 


T(0) = 4p?A3,/A? (38) 
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Here, as in Section 2.1, A is the Htickel determinant for the isolated molecule, from which in 
Aap its A-th row and B-th column are deleted and f is a measure of the interaction between the 
contact atom and its molecular neighbour. Again, what is the role of the LRF in this road to 
calculating/understanding molecular conductivity? 


4.2. The atom-atom polarizability as a conductivity indicator 


Going back to expression (7) for the diagonal elements of Coulson’s atom-atom polarizability, 
a general element mtag of this forerunner of the LRF can be written as 


feo 


ras = (1/7) | (an(iy)/AGy) Pay (39) 


—oco 


(the contour integral in the complex plane in Coulson’s formalism [31] is hereby reduced to an 
integral along the imaginary axis). The integrand of 714g thus turns out to be related to the 
transition probability at the Fermi level T, T4g(0) when the contacts are placed at atoms A and 
B of the molecule. An intimate relation between mag and T4,(0) thus exists; the precise 
connection was evaluated by explicit calculation of 774g and T,4g(0), some results being visual- 
ized in Figure 9. There, we depict both the 714g and T,4,(0) values for some linear acenes 
(benzene, naphthalene, anthracene and tetracene), taking always one atom as the reference 
atom. For the atom-atom polarizability, the reference atom is denoted by a green circle with its 
area proportional to the self-polarizability of that atom which is as pointed out in Section 2.1 as 
always negative. Black and red circles on the other atoms denote the atom-atom polarizability 
values corresponding to a perturbation on the reference atom (varying a,) on the (charge of 
the) considered atom B(q,). Black circles correspond to mag > 0, red circles to mag < 0. For 
Tag(0), a similar approach is followed. The empty green circle denotes the position of the 
first contact ‘A’ (zero ‘ipso’ transmission [71]), the magnitude of the black circles on the other 
atoms ‘B’ denotes the magnitude of the transmission when the second contact is placed at 
position B; note that no red circles arise as T4g(0), the transmission probability, always lies 
between 0 and 1. 


Figure 9 shows that the pattern in the two plots is completely analogous and leads to the 
conjecture that a positive atom-atom polarizability seems to be a necessary condition in these 
Kekulean benzenoids in order to have transmission for a certain configuration of the contacts 
on the molecule. If the areas of the circles are considered, no exact proportionality between mag 
and T4g(0) can be inferred, but they are correlated. 


This issue was further investigated in the next linear acene and pentacene (Figure 10), by 
varying the position of the first contact. For a fixed first contact, the highest transmission 
occurs when the second contact is at the atom with the highest 7t4g value involving the first 
atom. From that atom on, a monotonously decreasing probability is noticed in either direc- 
tions. The sharpest decline in T is in the direction of the neighbour that exhibits the lowest 
atom-atom polarizability. On the basis of these results, it can indeed be conjectured that a 
positive atom-atom polarizability is a necessary condition for transmission and that the ten- 
dencies between the two series of values are similar. 
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Figure 9. The atom-atom polarizability (left) and transmission probability at the Fermi level (right) for a single reference 
atom for benzene, naphthalene, anthracene and tetracene (Reprinted from [71] with permission of AIP Publishing). 


Figure 10. The atom-atom polarizability (left) and transmission probability at the Fermi level (right) for pentacene for 
variable reference atom (Reprinted from [71] with permission of AIP Publishing). 


Further analysis of the behaviour of the (A,z/ A)* function in the complex plane was done in 
the case of alternant non-singular hydrocarbons (for details see [71]). Coulson’s and Longuet 
Higgins’ pairing theorem shows that if A and B are drawn from the same partite set (all 


‘starred’ or ‘unstarred’ atoms, respectively) A3.,(z)/A (z)* is odd and if A and B are taken from 
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a different set this function is even. It is then easily seen that the odd function has no real part 
along the imaginary axis, leading to a negative (A4g/A)’ value along the y axis and zero at the 
origin, yielding a negative atom-atom polarizability and zero transmission at the Fermi level 
(T(0)). On the other hand, when A and B are drawn from opposite sets, Ayg/A is an even 
function yielding real values along the y axis. A positive 774g value results with either insula- 
tion or transmission, depending on whether A4g(0) is zero or not. Note that the analysis for 
A = B leads to the same result as for A # B but belonging to the same set yielding a negative 
Taa Value. The following overall conclusion, formulated as a selection rule, can be drawn: 
negative atom-atom polarizabilities for non-singular alternant hydrocarbons correspond to 
devices with insulation at the Fermi level. Conduction at the Fermi level requires, but is not 
guaranteed by, a positive atom-atom polarizability. 


The aforementioned properties were used as guiding principle in our later studies towards a 
chemical interpretation of molecular electronic conductivity [77] leading to a simple, back- 
of-the-envelope determination of quantum interference [78], thus bridging the gap between 
chemical reactivity theory and molecular electronics. 


5. Conclusions 


The LRF and its congener, the softness kernel, are now in a stage where many of their 
mathematical and physical properties are well understood. The possibility to evaluate, repre- 
sent and interpret them puts them on equal footing for their use in conceptual DFT with their 
already more traditional second-order companions, the chemical hardness and the Fukui 
function. In view of the ‘chemistry’ contained in the LRF kernel as shown some years ago, it 
is not unexpected, but it still remains to be unravelled whether they are major players in very 
fundamental issues pertaining to the electronic structure of matter as in Kohn’s nearsighted- 
ness of electronic matter principle, as well as in more applied fields where they are shown to be 
of great use to explore chemical compound space (through the alchemical derivatives) and to 
predict/interpret molecular conductivity. 
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